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Abstract 

One of the most challenging issues in QCD is the investigation of spontaneous chiral- 
symmetry breaking, which is characterized by the non-vanishing chiral condensate when the 
bare fermion mass is zero. In standard methods, one has to perform expensive lattice simu- 
lations at multiple bare quark masses, and employ some modeled function to extrapolate the 
data to the chiral limit. This paper applies the probability distribution function method to 
computing the chiral condensate in lattice QCD with massless dynamical quarks, without any 
ambiguous mass extrapolation. The results for staggered quarks indicates that the method 
might be a more efficient alternative for investigating the spontaneous chiral-symmetry break- 
ing in lattice QCD. 

1 INTRODUCTION 

Chiral-symmetry breaking is an important feature of quantum field theory with fermions. In QCD 
with massless quarks at low temperature or density, the system is confined and chiral symmetry 
is spontaneously broken. At sufficiently high temperature or density, there is a chiral phase 
transition from the confinement phase to the quark-gluon plasma phase or color-superconducting 
phase, where chiral symmetry is restored. There have been many phenomenological or QCD- 
inspired model studies of spontaneous chiral-symmetry breaking and chiral phase transitions, but 
lattice gauge theory (LGT) [1] should give the most reliable information from the first principles. 
These issues[2, 3, 4, 5, 6, 7, 8, 9, 10] have been extensively investigated on the lattice for a long 
time, but there are still many unsolved problems. 

The main question we would address here is how to get quantitative information on spon- 
taneous chiral-symmetry breaking from lattice simulations with dynamical quarks. When the 
bare fermion mass m is zero, the action is invariant under the global chiral transformation of the 
fermion field ip — > expect • r^^)ip, with r the generator of the chiral symmetry group, but the 
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bilinear operator jpip changes. A non-vanishing vacuum expectation value {ijjip)\m=o is the signal 
for spontaneous chiral-symmetry breaking, which is also associated with dynamical generation of 
fcrmion mass. 

Suppose the quarks are degenerate. The chiral condensate per flavor is 

(W) = (TrA -1 )/ (iV c V), (1) 

where N c = 3 is the number of colors, V is the number of lattice sites and A is the fermonic 
matrix. The trace is taken in the color, spin and position space. The expectation value "< ... >" 
is computed with the integration measure associated with the partition function 

Z = J [dip] [dtp] [dU] exp(-S) = J [dU] exp (S g + In det A) , (2) 

where U is the gauge link variable, S g is the gluonic action, and A = m + iA is the fermionic 
matrix. 

A popular way to calculate the chiral condensate is to use the matrix inversion technique. 
Another popular way is to use 

1 l NcV/2 2m \ 

3 — ^ 

with Xj the j-th positive eigenvalue of A. The disadvantage is that it requires the calculation of all 
eigenvalues of the Dirac operator. When the lattice volume is large, the computational task is huge 
and not so feasible. It is also very expensive to generate configurations with dynamical quarks at 
multiple fermion masses. To get the chiral condensate in the chiral limit lim m ^o h m v->oo(V'V ; }> 
it requires a m extrapolation, since < ipip > \ m =o = on a finite lattice. Such an extrapolation 
of (tpip) to the massless limit uses some modeled fitting functions (e.g., linear function, plus 
quadratic or logarithmic corrections). Unfortunately, such a process might not be well justified, 
and sometimes it gives arbitrary results, in particular when the system is at a phase transition 
point. 

In Ref.fll], an alternative, named the probability distribution function (p.d.f.) method of 
the chiral condensate was proposed to investigate spontaneous chiral-symmetry breaking in LGT 
with fermions. The p.d.f. method has been tested in the Schwinger modelfll], and applied to the 
investigation of the spontaneous P and CT symmetry breaking[12], theta- vacuum like systems[13] 
and the phase transition of SU(2) LGT at finite density[14], as well as SU(3) LGT with quenched 
staggered quarks[15]. 

In this article, we will compute the chiral condensate of lattice QCD with dynamical staggered 
quarks, using the p.d.f. method. In Sec. 2, we review the basic ideas. In Sec. 3, we present the 
numerical results. Some discussions are made in Sec. 4. 

2 FORMULATION 

Let us characterize each vacuum state by a and the chiral condensate by {ipi/})a- The p.d.f. of 
the chiral condensate in the Gibbs state is defined by [11, 15] 

P(c)=J2w a 6(c-(W) a ), (4) 
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with w a the weight to get the vacuum state a. P(c) tells us the probability to get the value c 
for the chiral condensate from a randomly chosen vacuum state. If the ground state is invariant 
under the global chiral transformation, then P(c) = 5(c). If the ground state breaks the chiral 
symmetry, P(c) will be a more complex function, depending on the group of global chiral symme- 
try. Therefore, from the shape of the function P(c) computed in the configurations generated by 
a chiral symmetric action with exact m = 0, one can qualitatively judge whether chiral symmetry 
is spontaneously broken. 

In quantum field theory with fermions, chiral-symmetry breaking is dominated by the prop- 
erties of the fermion fields under global chiral transformation. From Eq. (4), one can define [15] 
the p.d.f. of the chiral condensate for a single gauge configuration U: 



[d^]exp(-S f )S 

Pu(c) = 



J2 X <p( x )*i>( x ) 

N C V 



/[#][#] exp(-5/) 



(5) 



where in "J2 X i^(x)^(x)" , summation over color, spin and position indices is implied. 

For gauge theory with fermions, if the action is invariant under £7(1) chiral transformation, 
we show [15] from p.d.f. that one can use 

C (j ) - M U)) - 0(1) (6) 

to compute the chiral condensate in the exact chiral limit. cq(U) is the amplitude of the chiral 
condensate for configuration U in the m = case. z(j) is the j— th zero of Jo, with Jo the 
zeroth order Bessel function of the first kind. The average is taken over gauge configurations with 
dynamical fermions at m = 0. In the chiral-symmetry breaking phase, a plateau for C(j) = const. 
will develop, from which the chiral condensate in the chiral limit can be extracted. 

The relation between eigenmodes and chiral-symmetry breaking is clear: if chiral symmetry 
is spontaneously broken, i.e., cq(U) ^ 0, according to Eq. (6), Xj should scale as z(j)/V. In 
the infinite volume limit V — > oo, the eigenvalues relevant for chiral-symmetry breaking are those 
going to zero as 1/V, which is consistent with Banks-Casher. 

The advantage of Eq. (6) is that to extract the value of C (j) from a plateau, only a few smallest 
eigenvalues for configurations with dynamical quarks atm = are needed for this calculation. Of 
course, finite size analysis limy^oo C(j) remains to be done, as in all approaches. 



3 RESULTS FOR QCD WITH DYNAMICAL QUARKS 

In Ref. [15], we have displayed how the method works in quenched lattice QCD. Here we would 
like to present the first data for SU(3) LGT with dynamical staggered (Kogut-Susskind) quarks. 
The action is S = S g + Sf. 

ly c p 

S f = ^(x)A x ^{y), 
x,y 
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U^x)U u (x + fi)Ul(x + iy)Ul(x), 



A 



(7) 



where (5 = 2N c /g 2 , with N c = 3. In the chiral limit m = 0, there exists a U(l) subgroup of the 
continuous chiral symmetry in the fermionic action and Eq. (6) applies. All simulations are done 
on the V = 10 4 lattice using the method described in [16]. 

Figure 1 shows the data for (tpt/j) for flavor number Nf = 1 at (3 = 5.1069 obtained from Eq. 
(3) at 14 nonzero m. For conventional simulation methods, it is very expensive to do computations 
at so many quark masses. A linear function is also used to extrapolate the data of {ipi/)) atm/0 
to the chiral limit. At this coupling, the linear fit works well for wide range of m € [0.004,0.1], 
so that the result can be used for comparison. 

Figure 2 shows the result of C(j) for the same f3 and Nf using Eq. (6) and from only one 
simulation [16] at m = 0. There is a nice plateau for j € [5,40]. The result is consistent with that 
from Eq. (3) extrapolated to the chiral limit and Eq. (6). 

Figures 3-8 show respectively the results for N f = 2 at /3 = 5.0238, N f = 3 at /? = 4.9411, 
and Nf = 4 at (3 = 4.8590. In the chiral limit m = 0, the data for different flavors have the 
same plaquette energy E p = 1.2921, but different /?, due to the effects of dynamical quarks. Our 
results indicate that the p.d.f. method gives estimate for the chiral condensate, consistent with 
the conventional one. 



To summarize, we have extended the p.d.f. method to lattice QCD with dynamical staggered 
quarks for computing the chiral condensate in the exact chiral limit. The first results are very 
encouraging. This might be an alternative efficient method for obtaining quantitative informa- 
tion on the spontaneous chiral-symmetry breaking in QCD. There are several advantages: only 
calculations of a small set of eigenvalues of the massless Dirac operator are necessary; only one 
simulation at m = has to been done at one f}\ there is no need for m or A extrapolation. 

In some cases, i.e., when the system is at the criticality, fitting the data obtained from obtained 
from Eq. (3) with some modeled function might give arbitary results, while using Eq. (6) might 
hopefully give more reliable estimate. 
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Figure 1: (ipip)/4 as a function of m for Nf = 1 and (3 = 5.1069 using Eq. (3). 
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Figure 2: C(j)/4 as a function of j for Nf = 1 and /? = 5.1069 using Eq. (6). 
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Figure 3: (tpip)/4 as a function of m for Nf = 2 and f3 = 5.0238 using Eq. (3). 
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Figure 4: C(j)/A as a function of j for Nf = 2 and f3 = 5.0238 using Eq. (6). 
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Figure 5: (V>^)/4 as a function of m for Nf = 3 and = 4.9411 using Eq. (3). 
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Figure 6: C(j)/4 as a function of j for Nf = 3 and = 4.9411 using Eq. (6). 
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Figure 7: (Y>-0)/4 as a function of m for Nf = 4 and (3 = 4.8590 using Eq. (3). 
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Figure 8: C(j)/4 as a function of j for Nf = 4 and /? = 4.8590 using Eq. (6). 



9 



